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Abstract 

This paper discusses wetting and capillary condensation transitions 
on a line and a rectangular array of cylinders using an interface poten- 
tial formalism. For a line of cylinders, there is a capillary condensation 
transition followed by complete wetting if the cylinders are sufficiently 
close together. Both transitions disappear as the cylinder separation is 
increased. The dependence of the wetting phase diagram of a rectangu- 
lar array of cylinders is discussed as a function of the chemical potential, 
substrate-fluid interaction strength and surface tension. 

PACS: 68.45.Gd and 47.55.Mh 
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1 Introduction 



The wetting of planar surfaces is now well understood jj], ^, |§j . However, 
in many realistic situations, substrates are far from planar. A particularly 
important example is provided by porous media,whose wetting properties 
have implications for fluid flow [W, oil recovery |5| and also the probing of 
the fractal geometry of surfaces |6| . Our aim in this paper is to describe 
wetting and capillary condensation on lines and arrays of cylinders as a 
step towards understanding the properties of binary fluids in complicated 
geometries. 

For such complicated geometries rigorous theoretical methods quickly 
become intractable and approximations must be made. A profitable ap- 
proach has been to use an interface potential which replaces the density 
profile at a fluid-fluid interface by a sharp kink and uses a local surface 
tension jl], |^| . This approach is valid far from the bulk critical point and 
for wetting layers thicker than a few intermolecular spacings. It has the 
advantage that it is easily implemented, yet provides qualitatively correct 
phase diagrams. 

Cheng and Cole fjj and Napiorkowski et al. Q applied the interface 
potential approach to wetting in a corner, Darbellay and Yeomans || to 
wetting in a slit, and Robbins et al. to wetting on a line of slits. Dobbs 
et al. fill] extended the approach to treat two spheres, and subsequently 
a square array of cylinders jl^] in both the grand canonical and the more 
physically realistic canonical ensembles. In each case, sensible qualitative 
results were obtained for the phase behaviour although, as one might 
expect, subtle details of the interface position are not given correctly H. 

In Section 2 we use such an interface potential approach to study 
wetting on a line of cylinders. We consider the case of Van der Waals 
interparticle interactions. If the cylinders are close enough together for a 
capillary condensation transition to occur, the system undergoes complete 
wetting as A/i — » 0. Otherwise, the cylinders behave individually rather 
than collectively and the wetting transition is suppressed. 

In Section 3, a rectangular array of cylinders is considered. The phase 
diagram is determined as a function of the aspect ratio of the array, the 
chemical potential, the strength of the Van der Waals interactions, and 
the surface tension. Limiting cases in which this substrate reverts to the 
line of cylinders described in Section 2 and the square array considered 
by Dobbs and Yeomans [Jl2| are discussed. 

Although using the interface potential approach significantly simplifies 
the problem, it is still not easy to study the more complicated geometries 
which will model a porous medium more closely. For example, even for 
an array of spherical substrates, the non-linear differential equations that 
must be solved become two- rather than one-dimensional. Therefore, in 
Section 4, we describe a simpler way of modelling the different phases, 
which is amenable to extension to more complicated regular and random 
geometries. A comparison to the solutions of the interface potential ap- 
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proach allows us to assess its usefulness. 
Our results are summarised in Section 5. 



2 An Infinite Line of Cylinders 

Firstly, we consider wetting on an infinite line of identical cylinders ly- 
ing along the a;-axis with their axes parallel to the z-axis, as shown in 
Figure f . The cylinders have radius ro and their separation is L' . The 
relevant part of the grand potential per cylinder, per unit length in the 
z-direction, is taken to be 



4> = 4^ + ' jcrVI + W) + Am (l(x) - gf) | dx + W[l(x)]j 

(1) 

where l(x) is the interface position and the subscript x denotes differenti- 
ation with respect to x. 

The first term in equation (|l|) is the free energy of the liquid-gas 
interface, the surface tension a multiplied by the surface area. The second 
is a bulk term due to the excess cost of the adsorbed, unfavourable liquid. 
If pi and p g are the liquid and gas number densities, respectively, then 
the free energy per unit volume of the liquid phase over and above that 
of the gas phase is 

Am = (/Xc - p){pi - p g ) (2) 

with p* the chemical potential of the fluid in the system and p c the 
chemical potential at bulk liquid-gas coexistence. 

The final term models the interparticle interactions which, for non- 
retarded Van-der-Waals forces, can be written 

W[l(x)) = / / U\x',y(x'),r ,L']dydx' (3) 

JO Jl(x) 

with a disjoining pressure 

n(x,y,r ,L')=yf — ^ dr' (4) 



/ cylinder i |H 

where the summation is taken over all cylinders i. Note that the integrals 
in equations (|^) and (Q) are over the gas and substrate; all other inter- 
actions are either independent of l(x) or can be reformulated as integrals 
over these regions. The strength of the interactions is Wq — A/tv 2 where 
A is the conventionally defined Hamaker constant. 

The integral in equation (Q) cannot be performed analytically, in con- 
trast to the cases of a spherical substrate Jf3| and a cylindrical pore Jf^] , 
but the numerical result is well fit by a function 
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n [x ,y,r ,L<] « rk^F + ~w) (5) 

where U is the distance from the centre of cylinder i to the point r. The 
first term in the expression ^ is accurate at small distances from a cylin- 
der, when the substrate acts like a flat plane. The second gives the correct 
behaviour in the long distance limit. The fit (^) agrees with numerical 
integration of the disjoining pressure to within 10%, with the largest dis- 
crepancy occuring in the crossover region, at a distance of about ro from 
the surface of a cylinder. This discrepancy is unimportant because the 
contribution of the interactions to the total free energy is small compared 
to that of the surface energy at this distance. 

The grand free energy, 4>, can be minimised with respect to l(x) using 
the Euler-Lagrange formula, yielding a non-linear differential equation 

?k=)-Aii + Il[x,l(x),ro,L f ]=0 (6) 

1 ~t~ lx) J 

A solution where the liquid forms bridges between the cylinders (see Fig- 
ure 1) may be found if this is solved with boundary conditions l x = at 
x = and x = L'/2. 

To find the unbridged solution, where the interface wraps around each 
individual cylinder, it is necessary to formulate the problem in polar co- 
ordinates with the origin at the centre of a cylinder. The Euler-Lagrange 
minimisation then gives 




d_ ( oh \_^ m+Jl \ e m L n l{e) ^ n (7) 

The boundary conditions are lg = at 8 — and 8 = 7r/2, where lg = 
dl(8)/dd. 

The solutions to the differential equations (^|) and (^) can be found nu- 
merically using a relaxation method for different values of the parameters 
fi — Afiro/cy, L = L' /ro and Wo = Wq/ '(org). Once the interface profiles 
are known, the grand free energy of each phase can be calculated from 
equation (gj) using numerical integration techniques, allowing comparison 
of the free energies and determination of the stable configuration. 

The resulting phase diagram is shown in Figure 2. For L < 2 the 
cylinders are overlapping and the problem is not defined. The unbridged 
phase is stable at large /i, as expected. For 2 < L < L c (Wo), as fi is 
decreased there is a first-order phase transition to the bridged phase and 
then, as fi — > the interface unbinds to infinity, an example of a complete 
wetting transition. As the surface unbinds to large distances, it becomes 
flat and it follows from (0) and (@) that 
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~ 7T 2 rg Wl, , s . 
M ~ 2 1// 4 aS °° ® 
This shows that the line of cylinders is behaving, as expected, like a plate 
of effective thickness tttq/L' . 

For L > L c , capillary condensation does not occur as fj, decreases and 
the complete wetting transition is suppressed by the substrate geometry. 
For vanishingiy small Wo at fi = 0, the surface area is the only relevant 
quantity. The liquid-gas surface area per unit length for a single cylinder 
in the unbridged phase is 2nro while that of a bridged film is 2L c ro, giving 
L c {0) = 7r in agreement with the numerical solution. 

As Wo is increased at fixed fj,, the interfaces lie further from the sub- 
strates. For an approximately flat, bridged interface the surface energy 
is virtually unchanged by this. However, for the unbridged solution the 
surface area and hence the surface energy must increase as the interface 
moves. Thus the unbridged solution becomes less favourable for a given 
fx, as seen in Figure 2. 

3 An Infinite Rectangular Array of Cylin- 
ders 

An infinite number of lines of cylinders can be brought together to produce 
a rectangular array with inter-line distance D' , as shown in Figure 3. As 
D — D' /ro is reduced, this system shows a crossover from the behaviour 
of a line of cylinders to that reported in Dobbs and Yeomans jl^] for a 
square array. 

Two new phases might be expected to exist. The first of these consists 
of bridging between lines of cylinders as well as between the cylinders in 
one line. To find such a doubly-bridged solution it is necessary to use a 
polar coordinate system centred on an interstitial site such as point A in 
Figure 3. The Euler-Lagrange minimisation then gives 

= (v#l) + *' w - n l9 ''"" %r °- L ' ] " 9K ^?I = ° <9) 

where I' and 0' are the distance and angle to point (I, 9) from the centre of 
one of the four nearest cylinders. The boundary conditions are that lg — 
at 8 — and 8 — ir/2. Equation Q is solved numerically, as before. 

There is also a phase where the space between the cylinders is com- 
pletely filled with liquid. The free energy of this per cylinder, per unit 
length in the z-direction is 

$ full = Ap(L'D' - nr 2 ) (10) 

From comparisons of the free energies, the stable configuration can 
be found for given values of the four parameters L,D,/j, — Afj.ro/ a and 
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Wo — Wo/err^. Typical values of the Hamaker constant for a condensed 
phase are of order (0.4 — 4) x 10~ 19 J, while typical interfacial energies lie in 
the range (2 — 8) x 1CP 2 J m~ 2 |D| . If we suppose that the cylinders are of 
the same size as the grains in a porous rock, then ro ~ (1 — 100) x 10 _6 m, 
giving Wo ~ 10" 6 . 

Phase diagrams are plotted in Figures 4 and 5 and are discussed below. 

• Small L 

Figure 4 shows a cross-section through the phase diagram for L = 
2.2 and Wo = 2.5 x 10 -6 . For any finite D, the behaviour is no 
longer that of many separate horizontal lines of cylinders. Complete 
wetting at fi = is replaced by a transition to a full phase at /i > 0. 
For large D, this transition lies along the line D — 2/fi. As D 
decreases, the phase which is singly-bridged along the ^-direction 
becomes unstable, as expected, and the doubly-bridged solution is 
favoured. For D = L — 2.2, symmetry demands that singly-bridged 
solutions cannot be stable and we recover, as expected, a phase 
diagram topologically similar to that of Dobbs and Yeomans for a 
square array 113] . 

For D < L, as n is decreased the system jumps from being unbridged 
to singly-bridged along y to doubly-bridged to full, the order of in- 
creasing liquid volume. 

• Increasing L 

For L — 2.6, Wo = 2.5 x 10 -6 , the phase diagram is that shown in 
Figure 5. As L is increased, the quadruple point B in Figure 4 moves 
to lower values of \i until it coincides with point A, when the doubly- 
bridged phase disappears from the phase diagram. Moreover, the 
phase which is singly-bridged along x is stable only at increasingly 
high values of D as L is increased further until L = L*, when it 
becomes thermodynamically unstable for all D. L* is, as expected, 
approximately equal to L c , the critical value above which a line of 
cylinders does not undergo complete wetting, with small corrections 
due to the influence of the other cylinders in the array. 

• Increasing Wo 

When Wo is increased, at a given value of jj, the free energy of phases 
with more interface close to a surface will increase relatively more 
than those with less. Hence the transition lines move to greater 
values of n. For Wo = 2.5 x 10 -5 and L = 2.2, the phase diagram is 
like that in Figure 4 but point B moves to /i ~ 6.67 while point A 
moves to values of fi and D greater by ~ 3%. 

4 Simple Model 

The interface potential approach already involves considerable approx- 
imation. However, it still relies on the high symmetry of the system 
considered to produce tractable, one-dimensional, non-linear differential 
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equations. To better model porous media it would be desirable to be able 
to treat more complicated substrate geometries. Thus, we now consider 
a much simpler way of modelling the phases in the cylindrical array and 
compare the resulting phase diagrams with those obtained from the inter- 
face potential approach. We find only small discrepancies for physically 
realistic Wo, giving confidence that the simpler model will give qualita- 
tively correct phase diagrams for more realistic models of porous media. 

The approach is to approximate the interface shapes in the different 
thermodynamically stable phases by simple curves that can be handled 
analytically. The effect of the Van der Waals interaction is incorporated 
by assuming that where an interface wraps around a cylinder it lies a 
distance r from the centre of the cylinder, where r follows from the flat 
plane result 

r-r =?[m (11) 

V 6A m 

Where the liquid forms bridges between cylinders, the bridges are taken 
to have a radius of curvature a = a/A/j, = ro/fi, which follows from 
minimising the free energy of a bridge with respect to a. By assuming 
that the arcs of radius r and a meet tangentially, the interface shape is 
completely defined. 

The different phases can be modelled as shown in Figure 6. 

• Unbridged phase 

A circle of radius r, centred on a cylinder. 

• Singly-bridged phase 

An arc of radius r, centred on a cylinder, for 9 > Go, with 9 measured 
from the direction of bridging. An arc of radius a tangential to this 
at 6 = Oq , with -IP- = at x — L' 72, where cos Oq = , L , , . 

u ' ox i ' J 2(r+a) 

• Doubly-bridged phase 

An arc of radius r, centred on a cylinder, for Oq < 9 < 0\. Two 
arcs of radius a tangential to this, the first at 9 — 9o with = 
at x = L' 12, where cos#n = — r and the second at 9 = 9i with 
§ = at y = D'/2, where cos 9 1 = 

The grand free energy is taken to be 

$ = a x interface area + A^t x fluid volume (12) 

The stable state can be determined by comparing the free energy of the 
different phases. Despite the simplicity of the approach the phase dia- 
grams in Figures 4 and 5 are reproduced to within 1%. 

Note that a term representing the interaction energy has not been 
included explicitly in equation (J12|) . It can be estimated by, for example, 
considering only those parts of an interface which adhere to the substrate 
as making a contribution. However, for the unbridged solution the ratio of 
the surface free energy to the interaction or volume terms ~ \^Wo- For a 
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bridged phase the volume component of the free energy increases, but the 
contribution from the interactions remains of the same order of magnitude. 
Therefore, for the values of Wo considered here, the contribution of the 
interactions to the total free energy ~ 1% and can reasonably be neglected 
given the level of approximation already inherent in the approach. 

5 Discussion 

In this paper we have described wetting on a line and rectangular array 
of cylinders. For a line of cylinders a capillary condensation or bridging 
transition is followed by complete wetting as /i — » 0. If the cylinders are 
sufficiently far apart bridging does not occur and the wetting transition 
is suppressed: the cylinders are now behaving individually rather than as 
an effectively planar substrate. 

For an array of cylinders we have calculated the phase diagram as 
a function of the aspect ratio and reduced chemical potential. Several 
different capillary condensation transitions occur: to states bridged in the 
x or y directions; to a doubly-bridged phase or to a phase where the liquid 
completely fills the volume between the cylinders. 

The thin-thick transitions which correspond to wetting on a cylindrical 
substrate [H| are not included in this model where we consider an 
effective interface potential with a single minimum. The results of Dobbs 
and Yeomans jl^ for adjacent spheres indicate that including retarded 
Van der Waals terms in the potential to model such transitions would not 
substantially affect capillary condensation, while the thin-thick transition 
lines would essentially follow those for an individual cylinder. 

The phase diagrams have been obtained using an interfacial potential 
approximation which treats the interface as a sharp delineation between 
the two phases. The interface position is obtained as the solution of a 
non-linear differential equation. The feasibility of the approach is depen- 
dent on the high symmetry of the system we have considered. To better 
model porous media it is important to be able to treat more complicated 
substrates. Therefore we have tested a simpler approach where the shape 
of the interface in each phase is fed in as an assumption. The results agree 
very well with those obtained by the interface potential method giving us 
confidence that the simple approach will give realistic results for substates 
approximating more closely those found in porous media. 
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There have been micromodel experiments undertaken on regular two- 
dimensional arrays jl8| where capillary condensation is seen to play a role 
in the flow of a binary liquid mixture through the model porous medium. 
Our results apply to static configurations of the liquid but may supply 
insight into the flow properties. 
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Figure Captions 

Figure 1 The arrangement of fluid around a line of cylinders for a bridged 
phase. 

Figure 2 Dependence of the bridging transition of a single line of cylin- 
ders on the reduced Van der Waals interaction Wo'- solid line, Wo = 
0; long-dashed line, Wo = 2.5 x 10~ 6 ; shorter-dashed line, Wo = 
2.5 x 1(T 4 . 

Figure 3 Part of a rectangular array of cylinders, showing the doubly- 
bridged phase. 

Figure 4 The phase diagram of an array of cylinders with L = 2.2 and 
Wo = 2.5 x 1(T 6 . 

Figure 5 The phase diagram of an array of cylinders with L = 2.6 and 
Wo = 2.5 x 1(T 6 . 

Figure 6 The approximate geometries used to model the interface po- 
sition in the approach described in §4 : (i) unbridged phase, (ii) 
bridged phase, (iii) doubly-bridged phase, a — r //i. 
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